Magnetic resonance based apparatus and method to analyze and to measure the bi-directional flow regime in a transport or a production conduit of complex fluids, in real time and real flow-rate

ABSTRACT

An apparatus and a method is provided based on the Magnetic Resonance techniques to analyze and measure a uni- and/or a bi-directional flow regime of multiphase fluids in a transport and production conduct, in real time and flow-rates, based on a magnetic resonance analytical module, two magnetic prepolarization modules of variable effective length and a control computer of data acquisition and transfer, all associated to each other.

FIELD OF THE INVENTION

The present invention is related to the direct determination in theproduction line, of the flow regime in a multiphase complex fluid, bythe selective measurement by Magnetic Resonance (MR) of the velocityprofile of each separate phase, in a transversal slice to the flowdirection. The present invention is specifically related to a method andan apparatus to measure, in real time and real production volumes, theindividual component flow-rates, the relative proportions of thedifferent phases that compose the complex fluid and the spatial profileof velocities in the fluid vein, in both uni and bidirectional flowingfluids. In particular, the present invention is applied preferably tothe analysis of the flow in oil industry production veins, as well as inother conducts that transport diverse types of fluids, including, forexample, the cases where water and petroleum move in oppositedirections. The analysis is performed on a non-invasive andnon-destructive way and independently from the state of the mixture,this means for example, if it is composed by petroleum, water and gas inseparated phases or in a phase in emulsion.

Other applications in the mining, food, chemical, pharmaceuticalindustry, etc., are also possible, being for instance one fluid composedof miscible and/or non-miscible pastes, or mixtures of liquids that havesolids in suspension.

BACKGROUND OF THE INVENTION

Many solutions exist in the literature using Magnetic Resonancetechniques to analyze directly on the liquid pipe and in real time,different properties of the transport and/or production of fluids. Wecan for instance mention:

1. Up from the measurement of the relaxation times T₁ and T₂ and fromthe coefficient of Molecular diffusion D, the petrophysical propertiesof the extraction fluid are measured in real time and in the realconditions of reservoir, including:

-   -   a. the level of contamination of the fluid due to, by example,        filtrations of water, mud, etc.;    -   b. viscosity and the relation water/gas/petroleum.

These measurements are made in surface as well as in the pit(continental or offshore) and eliminate ambiguities associated withsamples extraction as well as the corresponding transport procedures.Some analytical methods based on the Magnetic Resonance and itsassociated devices are described, for example, in the U.S. Pat. No.6,111,408, “Nuclear Magnetic Resonance Sensing Apparatus and forTechniques down hole Measurements”, by T. Blades et al.; U.S. 6.737.864B2, “Magnetic Resonance Fluid Analysis and Method”, by M. G. Prammer etal.; U.S. 6.825.657 B2, “Magnetic Resonance for Method CharacterizingFluid Samples Withdrawn from Subsubsurface Earth Formations”, by R. L.Kleinberg et al.; U.S. 6.891.369 B2, “Magnetic Resonance Method andLogging for Apparatus Fluid Analysis”, by M. D. Hurlimann; U.S.2005/0040822 A1 “Multi-measurements NMR Analysis based on MaximumEntropy”, by N. Heaton and 2006/0122779 U.S.A1 “Interpretation forMethods NMR Diffussion-T2 Maps” by Chang Cao Ming and N. Heaton and thetherein contained references

2. Another group of tools, also based on the Magnetic Resonancetechnique, has been divulgated, to determine the cut (or proportion ofpetroleum and water) and the flow in the fluid vein. On the one hand,the determination of petroleum and water cut is generally is carried onby the Magnetic Resonance signal weighted by the individual spin-latticerelaxation times T₁ of the fluid. In particular, for the case of amixture of petroleum and water, water T1 differs from those ofpetroleum. Moreover, it is also possible to measure the cut of the lightand heavy elements of petroleum, since its T₁ values differ enough toisolate the different MR signals. See by example the patent No. U.S.4.785.245. On the other hand, for the flow-rate measurement, threedifferent basic principles can be grouped:

The measurement of the fluid flow-rate through the measurement of the“flight time” of the fluid between two Magnetic Resonance spectrometers:(or between two sensors of a same spectrometer). See the U.S. Pat. No.6,046,587 “Measurements of Flow Fractions, Flow Velocities and FlowRates of a Multiphase Fluid using NMR Sensing” and the U.S. Pat. No.6,268,727 “Measurements of Flow Fractions, Flow Velocities and FlowRates of a Multiphase Fluid using ESR Sensing”, by J. D. King, Q. Norand A. De los Santos, which disclose a sensor that uses at least two MRspectrometers or one MR and another Electron Paramagnetic Resonance one.The basic principle of the measurement methodology is based on what isknown as the “flight or passage time” of MR-excited fluid nuclei betweenboth MR spectrometers. Another variant is the US Patent Application No.2004/001532, by M. Ramia, D. J. Pusiol, C. A. Martin, E. Fried and R.Garnero, “Method and procedure to measure fluid flow and fluid fraction,and equipment used to that end”. In this case there is only oneelectronic part, shared by two sensorial coils, being the operationprinciple the same as the one described before; namely, the speed of theflow through flight time of the water and petroleum molecules in thespace between both sensorial coils is separately measured. Those arelittle practical and expensive to be implemented, as well as ofdifficult application in petroleum fields with rigorous climaticconditions generally.

Another methodology is based on the spatial encoding of the flowvelocity by a magnetic field gradient in the direction of the flow. Thegradients (static and/or electronically pulsed) are employed to modulatethe protons spin precession phase. It means that the spatialcodification is made in what it is known as the Laboratory Frame. Theflow meter with fluid phase separation that uses the pulsedelectromagnetic field gradient is disclosed, for instance, by the U.S.Pat. No. 6,452,390, by E. Wollin, “Magnetic Resonance Analyzing FlowMeter and Flow Measuring Method”. This method has the disadvantage that,at the speeds at which the protons commonly move in the magnetic fieldwhen the fluid flow-rate is measured under reservoir conditions, is toofast. Therefore, the application of fast rise time magnetic fieldgradients switching on and off is of difficult technologicalimplementation. This means that this methodology is generally restrictedto measurements at relatively low flow-rates. Another variant whichavoids high electric currents by including a permanent longitudinalgradient field is described in the US Patent Application No. US2006/0020403, by D. J. Pusiol “Device and Method for real time directmeasurement of the Flow-Rate of a Multi-Component Complex Fluid”. Thisinvent divulgates a flow meter and a cut sensor constituted by only onecoil associated to a magnet of slightly oblique flat polar faces. Thedevice generates a constant magnetic field, in addition with a magneticfield gradient in the direction of the flow. The spatial codification ofthe temporal position of the resonant nuclei is done by a linearmagnetic field gradient in the volume that occupies theexcitation/detection Magnetic Resonance coil. For high flow speeds, thisgradient must be increased to reach the necessary effectiveness in thespace codification process of the protons that compose the circulatingcomplex fluid. But the increment of the magnetic field gradient strengthimplies a broadening in the MR line width. Therefore, the MR signal tonoise ratio deteriorates. The maximum limit of the flow that is possibleto be measured by the mentioned invention is determined by the maximumvalue of the magnetic field gradient that is possible to apply to thefluid, before the signal is deteriorated so that the required precisionis not anymore reached for this particular use.

All the above described solutions measure only the average of flow-ratesof the fluid components, but to evaluate the efficiency of the methodsof extraction, pumping, transport, water injection in secondaryextraction, etc., it is necessary to be able to evaluate, in the lineand in reservoir conditions, the flow regime of the fluid. Therefore, itis also necessary to know the velocity profile of each individualcomponent of the fluid in spatially selected volume elements.

SUMMARY

The above drawbacks and difficulties are solved by this invention, asthe invented method and device allow real time measurement of theprofile of the flow-rate and proportions of the different componentsconforming a multicomponent complex fluid in a non-invasive andnon-destructive way, i.e. for example, regardless of the fact of itbeing formed by oil and water in separate phases or as an emulsion.

It is to borne in mind that this invention can also be applied to othertechnological fields such as, for example, that of the dairy industry,for determination of the milk-yogurt ratio, food and extruded organicmaterials, etc.

More specifically, this invention is related to a method and device fora Magnetic Resonance based apparatus and method to analyze and tomeasure a bi-directional flow regime in a transport or production pipeof complex fluids, in real time and flow.

The analytical magnetic resonance based method to be applied in the oilindustry is the following one:

The circulating fluid comes into a region of pre-polarization to codifythe component or phase whose relative proportion in the fluid andvelocity profile is to be measured or analyzed. For that purpose we usethe properties of longitudinal relaxation—characterized by T1—of eachcomponent of the fluid, the media speed of the fluid and the length ofsuch magnetic field.

Secondly, the properly polarized spins come into a region in which atransversal (and if not, a longitudinal) magnetic field has beengenerated, which also has a field gradient—properly designed accordingto the symmetry of the velocity profile to be measured—that can bereoriented in order to obtain different projections of the spins densityof the previously selected components and their relative speed.

Thirdly, an adequate excitation/detection pulse sequence measures andcharacterizes the flow regime by the following sequence:

The spatial selection of a volume element located inside the excitingcoil is made by a resonant radiofrequency soft pulse.

Once the volume element to be measured is selected, the average speed ofthe previously selected component in the prepolarization step, containedin such selected volume element, is measured. The slope of the echoeslinear decay characterizes the volume element velocity of the selectedfluid component. Such lineal decay in the spin echo signal of both thevolume and the selected fluid component refers essentially to theportion of said volume element that stays in the volume of theexciting-receiving coil. In the temporary ends of said decay dominatethe relaxation processes, and as such, this is characterized by anexponential component in the temporal evolution of the spin echointensity.

For each excitation frequency, or volume of selected fluid, steps a andb are repeated till obtaining a matrix whose elements represent thespeed and respective proportions in different spatial positions,preferably in the transversal slice of the fluid direction.

Finally, the previous steps are repeated for:

Different magnetic field gradient directions and/or intensities (spatialcodification).

Different effective lengths of the prepolarization field.

Different polarization fields (right or left) in order to make abi-directional fluid flow characterization.

The apparatus invented here to analyze and to measure the uni andbidirectional flow regime is a measurement and analysis module ofmagnetic resonance, two magnetic prepolarization modules of variableeffective length and a control computer of acquisition and datatransfer, all associated to each other. The measurement and analysis MRmodule is disposed symmetrically left and right with the twoprepolarization modules. The fluid conduct passes through that modulesarray.

The magnetic resonance analytical module generates a homogeneous,thermally compensated and constant magnetic field, with addition of aspace location gradient, constituted by a plurality of magnetic bars ofdifferent sizes and materials, that constitute a plurality of discs ofvariable thicknesses, piled up and associated among themselves.

The protons prepolarization magnetic field modules of variable effectivelength consists of successive discrete Halbach-type cylindricalsegments, that can be piled and associated to each other to select thecomponents of the fluid with successive minor spin-lattice relaxationtimes.

BRIEF DESCRIPTION OF THE DRAWINGS

The invention shall be best understood when read in conjunction with thefollowing figures:

FIG. 1: Illustrates a simplified block diagram of the preferredembodiment of the present invention.

FIG. 2: Illustrates another preferred embodiment, as it has been formedwith two measurement and characterization of the velocity profiledevices in the extremities of a sole prepolarization region of variableeffective length.

FIGS. 3 a)-3 c): Illustrate the previous art for a magnet of the Halbachtype with disc form, constructed with properly located discretecross-section rectangular magnetic bars.

FIG. 4: a) Illustrates a way to arrange the bars as described in FIG. 3;b) illustrates a convenient way to pile up those discs shown above.

FIGS. 5 a) and 5 b): Illustrate an embodiment of the previous art forthe homogenization of the internal magnetic field of the Halbachcylinder, which is commercially applied in Magnetic Resonance Imagingapplications.

FIG. 6: Illustrates the preferred embodiment for a measurement andcharacterization set of the velocity profile constituted by anadjustment of magnetic bars that compose the permanent magnet of thesensorial element, which has “Shimming” temperature compensation, andgradient generation in the x direction.

FIGS. 7 a)-7 c): Illustrate the preferred embodiment for the magnet ofvariable effective length for the prepolarization of the hydrogen nucleiof each component of the complex fluid.

FIG. 8 a) shows the succession of echoes that is obtained from theapplication of the CPMG sequence, for petroleum circulating at avelocity of 1.1 m/s; b) it shows the average fluid velocity of the −B/Aratio; where B is the slope and A is the y-axis interception at t=0 ofthe linear fit of CPMG echoes.

FIG. 9: Illustrates the MR measurement and analysis procedure todetermine the transversal velocity profile and the relative proportionsof the components of a complex fluid.

DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS

FIG. 1 shows a simplified block chart of the preferred embodiment ofthis invention. The complex fluid flows through the production vein 1 inany of the two possible directions; a phase can also flow in onedirection and the other/others in the opposite one 5. If the flow of thephase, whose regime needs to be measured or characterized, flows fromthe left to the right, for example, it comes first into theprepolarization magnet of effective variable length to polarize thespins of the phases that come from the left 3; then the fluid comes intothe measurement and characterization set of the velocity profile 2. Thefollowing step 4 produces no effects in this situation, as it isadjusted to an effective zero polarization length. The magnetic fieldgenerated in the first step 3 allows choosing the phase of the fluid tobe measured in the following step 2. This operation is made by therelation between the average velocity of such phase in the pipe, thevalue of the relaxation time T1, and the effective length that the fluidflows in such prepolarization magnetic field. For fluids with more thantwo phases, this procedure is repeated as many times as necessary, eachtime with different effective prepolarization lengths. We will furtherreturn to the detailed description of this point. Once the cut iscodified as explained before, the fluid comes in the set of localmeasurement of the velocity profile 2. In a preferred embodiment, thisstage consists of a magnet that generates a magnetic field in atransversal direction to the flow, able to be used in experiments ofmagnetic resonance, which have also a magnetic field gradient in aCartesian direction for a preferred application, or in a polar (orradial) direction for another preferred application. Such Cartesiandirection may be oriented by the rotation of the measurement device, orby the disposition of many electromagnets arranged in a convenient way.Inside the magnet are the radiofrequency coils used to excite the nucleiof the prepolarized phases in step 3, and to receive the magneticresonance signal. The sequence of events is the following:

First, all the nuclei that are in a desired portion of the total volumecontained in the radiofrequency coil are excited. This volume isdetermined by a first radiofrequency pulse, with pure spectralcharacteristic, or narrow band, also known as “soft pulse”.

Immediately, a sequence that reads the quantity of the previouslyselected spins still remaining in the volume of the detection coilstarts. Such reading is preferably made by a sequence of refocalizationpulses of the magnetization and also by the measurement of the generatedspin eco intensity. Said refocalization pulses may be narrow rangepulses (soft pulses), or wide range pulses (hard pulses). A preferredsequence of pulses is known as spatially selective CPMG (see forexample, R. Kimmich, “NMR—Tomography, Diffusometry, Relaxometry”,Springer-Verlag, Berlin—Heildelberg (1997)). As it later will bedescribed, in a range of intervals between the reading pulses, there isa lineal relation in the signal intensity decay of the spin eco. Thisinvention also proves that the slope of such lineal relation isproportional to the fluid local velocity. This process is repeated forseveral frequencies of the spatial selection pulse and severalorientations of the magnetic field gradient with respect to asingularity in the symmetry of the transversal velocity profile of thefluid

If necessary, other magnetic field gradient in the fluid direction canbe added in order to allow the selection a thick slice of the totaltransversal volume contained in the radiofrequency coil.

FIG. 2 shows another preferred embodiment. It is composed of twoanalysis modules of the flow regime (one for each direction of thefluid) located at both lateral ends; and in the center we reach a solesection of prepolarization. In this way, with two spectrometers we canmeasure simultaneously both fluid directions. Further we will describein detail each one of the parts and/or devices of the preferred processand the preferred set. To do that we revise a part of the previous art,especially that concerning magnets.

There exist in the previous art preferably cylindrical permanent magnetswith tubular geometry, based on polarized rare earth alloy segments. Themost known and commercially available are the so called “magicalcylinders”, standing out among them the Halbach magnets. The idea ofHalbach was to use polarized segments, to form the cylinder or ring, insuch a way that their magnetizations should fulfill the theoreticalrequirements for the angle (γ₁=(1+k_(i))β_(i)) (See, K. Halbach, “Designof permanent multipole magnets with oriented rare earth cobaltmaterials”. Nucl. Instr. Meth., 169, 1 (1980); and K. Halbach, “Strongrare earth cobalt quadrupoles”. IEEE Trans. Nucl. Sci. 26, 3882 (1079)).In the FIG. 3 a) a previous art magical Halbach cylinder composed by 8segments of trapezoidal shape 10, which generate in the inner volume asufficiently homogenous cross magnetic field 12 for uses in MagneticResonance experiments. This type of magnets have been improved,including a configuration to compensate thermal changes in the magneticfield (to see for example, M. G. Prammer and P. Masak, “Temperaturecompensated magnetic circuit”. Patent Application Publication U.S.2003/0210050 and Ms. G. Prammer and P. Masak, “Temperature compensatedmagnetic field apparatus for NMR measurements” U.S. Pat. No. 6,577,125B2). In a recent work H. Raich and B. Blümber, “Design and constructionof to Dipole Halbach Array with to Homogeneous Field from identical BarMagnets: NMR Mandhalas”. Concepts Magn. Reson. Part B (Magn Reson Eng),23B, 16 (2004) propose, as shown in FIG. 3 b), to use conventionalmagnetic bars 14, adequately arranged and oriented 16 in order to obtaina transversally oriented magnetic field. By this design, the authorsobtain a magnet of relatively small weight than that of a conventionalHalbach, although with smaller homogeneity in the useful space. FIG. 3c) illustrates the relative angles between the bars to obtain the mosthomogenous possible internal field, in accordance to the equation:

γ_(i)=(1+p)β_(i)

FIG. 4 illustrates a previous art embodiment proposed by Raich andBlümber. A plurality of rectangular magnetic bars 20 are arrangedtogether with two discs of non-magnetic material 22, fixed by the screws24 also of non-magnetic material. These discs 26 are piled up and heldby non-magnetic screws 28, as illustrated in FIG. 4 b), to form acylindrical permanent magnet. Before using Magnetic Resonanceapplications, it is necessary to homogenize and thermally compensate themagnetic field in the inner region of that cylindrical device.

FIG. 5 illustrates an embodiment of the previous art (see MAGNETICASOFT,Nice, France, “MRI Halbach 0.23 T magnet”, available inhttp://www.maxneticasoft.com/ of a magnet designed to be used in imagesdevices by Magnetic Resonance (MRI). The useful space for images byMagnetic Resonance 30 is surrounded by two concentric cylinders. Theinner cylinder 32 generates the main component of the field and theouter cylinder 34 is slightly rotated with respect to the first one andserves on one hand to increase the magnetic field in the zone 30, and onthe other hand to homogenize the magnetic field.

FIG. 6 illustrates a preferred embodiment for a set of measurement andcharacterization of the velocity profile. This device is constituted byan array of magnetic bars which compose the permanent magnet of themagnetic resonance analysis device. The rectangular bars array isdesigned to compensate thermal changes in the magnetization, innermagnetic field shimming, and the magnetic field gradient generation.Such bars array form discs of variable thickness with the followingproperties:

In the zone of interest 100 (delimited by the dotted circle in FIG. 6) amagnetic field Bx 111 is produced and a gradient of transversal magneticfield Gx 131 superposed to Bx, both having the following properties andcharacteristics:

Gx is transversal to the flow direction.

Bx is sufficiently homogenous and thermally stabilized for MagneticResonance experiments.

This field is generated by a plurality of magnetized bars 110, 120, 130and 140, which are properly arranged to produce respectively themagnetic field Bx 111, the homogenization effect or magnetic fieldshimming, the Gx Gradient generation 131 and the thermal compensationfield 121.

The shimming effect is double. On one side the bars 120 located in thediagonals contribute to compensate magnetic field in homogeneitiesproduced by the remnant magnetic quadrupole moment; on the other sidehigher order moments are compensated by a second set of bars arranged inthe directrix between the principal magnetic field bars array 110.

The geometrical arrangement and materials of magnetic bars 110 and 140allow the thermal compensation effect on Bx. It is well known that themagnetization produced by magnetic materials generally present anegative thermal coefficient. Moreover, the same effect is produced bythe thermal expansion of the bar array container. The basic principleconsists of the generation the Bx field by the superposition of twofields of opposite direction, Bx₁ and Bx₂, each one formed by aplurality of properly arranged magnetic bars. The thermal coefficientsof the magnetic bars with which are generated both fields have the samesign but substantially different values. Moreover, the magnetizationgenerated in the center of the magnetic set for Bx₁ is substantiallystronger than the one for Bx₂. For example, and without losinggenerality in the invention, the bars generating Bx₁ are chosen of analloy of Sm₂Co₁₇ (temperature Coefficient α=−0.03%/° C.) and the barsthat generate Bx₂ are chosen of Ferrite 8 (Coefficient of temperatureα=−0.20%/° C.). Furthermore, if the temperature of the magnetic setincreases, the diminution of the magnetization in the work zone 100 isproduced by i) the diminution of the magnetization of the Sm₂Co₁₇ bars,and ii) by the thermal expansion of the packing disks of the bars array.As the magnetization produced respectively by Sm₂Co₁₇ and Ferrite 8 areopposite, and its thermal coefficient presents the same sign, theindividual thermal changes in both magnetizations compensate the totalinternal magnetic field. Of course, in the calculation of the relativestrength of Bx₁ and Bx₂ the thermal expansion of the packing discs mustbe taken into account.

Bx can be transversal or longitudinal to the flow direction. Further,and without losing the generality of the invention, we will refer to Bxindicating the magnetic field transversal to the flow direction in thefluid vein.

The device, which generates the homogeneous and thermally compensatedmagnetic field, is obtained by piling up a plurality of discs asdescribed in FIG. 6, generally each one of different thickness. Inaddition, it is also possible to generate the spatially homogeneousmagnetic field for this invention by an arrangement of magnetic bars ofsuitable length. To the practical effects, to obtain the best possiblelongitudinal homogeneity, and without supposing a loss in the generalityof the invention, we refer further to the fact that the magnetic fieldof the measurement and characterization of the velocity profile device,is generated in the zone of interest 100 by a plurality of discs ofvariable thickness, piled up with variable separation between eachother.

All the magnetic system is surrounded by a material of high magneticpermeability 180. By this media an increment in the concentration ofmagnetic flow lines is produced in the zone of interest and furthermoreit contributes to obtain a better homogenization of the magnetic fieldin volume 100. This device contributes, in addition, to protect and toencapsulate the complete flow regime analytical system.

The system of magnetic bars used to generate the magnetic fields inapplication in Magnetic Resonance, arranged in discrete space positions,allows further constructing lighter embodiments than the conventionalpermanent magnets, allowing diminishing magnet construction costs, aswell as in the mechanical movements of the device parts. It is furtherpossible to interconnect electrical and/or electronic elements in theinterstitial spaces between the magnetic bars, diminishing in this waythe final volume and weight of the device.

Another embodiment consists in locating a plurality of magnetic bars130, generating the spatial localization gradient 131, in anothercylinder located in the external as well as internal part with regard tothe main magnetic cylinder. Both cylinders, the gradient one and theprincipal magnetic one, can rotate the one with regard to the other,producing the suitable reorientation of the Gx gradient 131.

The fluid vein is introduced in a measurement and characterization setof the velocity profile by a non-magnetic tube 170 of arbitrarygeometry, preferably of stainless steel, but with a radiofrequencytransparent section, preferably made of a fiber glass tube approved bythe oil industry, like the compound 316 SS. If the Bx field 111 istransversal, in this invention is preferably used a solereceiving/exciting coil of the solenoidal type 160. If another preferredembodiment is chosen, in which the magnetic field Bx 111 is longitudinalto the direction of the flow, the preferred receiving/exciting coil 160is of the “birdcage” type. A skilled technician can design and constructthese excitation and detection elements with the adapted bibliographicaid like for example: Jianming Jin, Electromagnetic Analysis and Designin Resonant Magnetic Imaging, pp. 193, CRC Press LLC (1999)

The grooved tube 150, constructed preferably in copper, contributes tothe electromagnetic shield of the radio frequency system.

A not in the illustrated mechanical system figures, is in charge ofrotating the set of measurement and characterization of the velocityprofile in all the positions necessary to obtain the complex matrix ofdata, which can be possibly represented graphically in image form.

A non illustrated Magnetic Resonance hardware and software, which iseither located in the proximity or away from the device of measurementand characterization of the velocity profile, but connected by signalreception and emission cables, is in charge of manipulating the radiofrequency signals and the pertinent digital information.

FIG. 7 illustrates a preferred embodiment for the magnet of variableeffective length for the prepolarization of the hydrogen nuclei of eachcomponent of the complex fluid. The device is constituted by discreteHallbach-type magnetic discs (see FIG. 4) of variable thickness andseparation 200, which can be rotated one with respect to the other by anot illustrated mechanical device. The effective polarization length ofthe magnet is determined by the relative number of discs whose internalfields are oriented in opposite directions. Obviously the maximumpolarization length is reached when the individual magnetic fields ofall the discs are oriented in the same direction, see FIG. 7 a). Thearray providing a prepolarization field of zero effective length is thatwhere the discs are oriented with equal number in opposite directions,see FIG. 7 b). Finally, FIG. 7 c) illustrates an arrangement of discsthat generates an effective polarization field of approximately onethird of the maximum field. A medium skilled technician can design andconstruct this set in accordance to the expected symmetry of the measurevelocity profile.

Measurements of average speed of the fluid and proportions of phasesMethods: the preferred method for measuring the average velocity of thetotal flow, as well as the proportions of phases present in a multiphasefluid with at least two phases that produce a detectable MR signal,consists of using a pulse sequence known in the previousstate-of-the-art like CPMG. The CPMG sequence is a sequence of pulses inthe form π/2−π− . . . −π, in such a way that the magnetization that willbe rotating on the x-y plane (detection plane) will be constructed withthose portions of the sample that were originally excited by the firstpulse of the sequence (the pulse of π/2) and that induces anelectromotive force in the detection coil. We define the time origin ofthe sequence (t=0) at the moment at which the pulse π/2 starts to beapplied. All the volume of the fluid inside the detection coil isexcited during such a short time interval.

In order to produce a better understanding of the invention, withoutsetting aside the general concept, the following analysis is restrictedto the simplest case in which the fluid is composed by two phases, forexample one of them is formation water and the other is petroleum oroil. In addition, in the present example, we will consider the case inwhich the production tube is completely filled with fluid during themeasurement process (at least, the total volume of the pipe must becompletely filled with fluid at the place where the MR sensor islocated).

So, in t=0 the total volume of the fluid inside the MR sensor will be:

V _(T)(0)=V _(A)(0)+V _(B)(0)  (1)

where V_(A) and V_(B) represent the volume occupied by each one of thetwo phases present in the fluid.

If the fluid moves on at an average velocity v, the volume of samplethat will produce a detectable MR signal at time t will be:

V _(T)(t)=V _(T)(0)−Svt  (2)

where S is the cross-section of the detection antenna.

-   -   At the time t, the following equations must be satisfied:

V _(A)(t)=V _(A)(0)−S _(A) vt=S _(A)(L−vt)  (3)

V _(B)(t)=V _(B)(0)−S _(B) vt=S _(B)(L−vt)  (4)

S _(A) +S _(B) =S,  (5)

where S_(A) and S_(B) represent the effective cross-sections inside thevolume of the detection coil, corresponding to the phases A and B,respectively, and L represents the length of the detection coil.

Dividing the equation (3) by (4) we obtain:

$\begin{matrix}{\frac{V_{A}(t)}{V_{B}(t)} = \frac{S_{A}}{S_{B}}} & (6)\end{matrix}$

The previous equation evaluated at time t=0 and at time t implies:

$\begin{matrix}{\frac{V_{A}(t)}{V_{B}(t)} = \frac{V_{A}(0)}{V_{B}(0)}} & (7)\end{matrix}$

This means that, in the present approach, at any moment of time withinthe interval of duration of the CPMG sequence, the relative proportionsof volumes of both phases are equal to they were at the moment t=0 whenthe pulse sequence began. So, the following expressions are easilydeduced from the equations (3), (4) and (7):

$\begin{matrix}{{V_{A}(t)} = {{V_{A}(0)}\left( {1 - {\frac{v}{L}t}} \right)}} & (8) \\{{V_{B}(t)} = {{V_{B}(0)}\left( {1 - {\frac{v}{L}t}} \right)}} & (9)\end{matrix}$

On the other hand, the intensity of the CPMG signal at time t will begiven by the following expression:

$\begin{matrix}{{{I\left( {v,t} \right)} = {{{I_{A}\left( {v,t} \right)}^{- \frac{t}{T_{2{eff}}^{A}}}} + {{I_{B}\left( {v,t} \right)}^{- \frac{t}{T_{2{eff}}^{B}}}}}}{{where}\text{:}}} & (10) \\{{I_{\alpha}\left( {v,t} \right)} = {{{M_{\alpha}(v)}{V_{\alpha}(t)}} = {{f_{\alpha}(v)}{M_{\alpha}(0)}{V_{\alpha}(0)}\left( {1 - {\frac{v}{L_{eff}}t}} \right)}}} & (11)\end{matrix}$

In equation (10), v represents the average velocity of the flow andT_(2eff) ^(α) represents the effective spin-spin relaxation time of thecomponent α. On the other hand, in equation (11), M_(α)(v) representsthe magnetization per volume unit of the component α circulating at anaverage velocity v, M_(α)(0) represents the magnetization per volumeunit of the component α at rest (v=0), V_(α)(0) is the volume of thecomponent α inside the MR sensor at the beginning of the CPMG sequence(t=0), and L_(eff) represents the effective length of theexcitation/detection coil. The coefficient f_(α)(v) is calledprepolarization factor and is given by:

${{f_{\alpha}(v)} = {\frac{M_{\alpha}(v)}{M_{\alpha}(0)} = \left( {1 - ^{- \frac{L_{pol}}{{vT}_{1}^{\alpha}}}} \right)}},$

where T₁ ^(α) represents the longitudinal or spin-lattice relaxationtime of the component α, and L_(pol) represents the length of theprepolarization stage, which means, the total length traveled by thefluid before entering to the MR sensor.

By the equation (11), the equation (10) can be written as follows:

$\begin{matrix}{{I\left( {v,t} \right)} = {\quad{\left\lbrack {{{M_{A}(0)}{f_{A}(0)}{V_{A}(0)}^{- \frac{t}{T_{2{eff}}^{A}}}} + {{M_{B}(0)}{f_{B}(v)}{V_{B}(0)}^{- \frac{t}{T_{2{eff}}^{B}}}}} \right\rbrack \left( {1 - {\frac{v}{L_{eff}}t}} \right)}}} & (12)\end{matrix}$

Of course, the expression (12) only has sense for those moments of thetime t that satisfy 0≦1(v/L) t≦1.

Determination of the average speed of the fluid: An importantobservation of the equation (12) is that, for values of t thatsimultaneously satisfy the conditions t/T^(A) _(2eff)<<1 and t/T^(B)_(2eff)<<1, this expression is reduced to:

$\begin{matrix}{{I\left( {v,t} \right)} = {\left\lbrack {{{f_{A}(v)}M_{0}^{A}{V_{A}(0)}} + {{f_{B}(v)}M_{0}^{B}{V_{B}(0)}}} \right\rbrack \left( {1 - {\frac{v}{L_{eff}}t}} \right)}} & (13)\end{matrix}$

The expression (13) shows that, if we make a fit of the amplitude of theechoes corresponding to the linear region of the CPMG sequence by anexpression such as I(v,t)=A+B t, the quotient −B/A, between the slopeand the y-axis intercept at t=0 of the fit, is directly proportional tothe average velocity of the fluid being the proportionality constantequal to the reciprocal of the effective length of theexcitation/detection antenna.

Determination of the fraction of components with prepolarization stageof fixed length (Method 1): The expression (13) shows that the y-axisintercept at t=0, obtained from the fit of the region where theintensity of the CPMG signal varies linearly for a given averagevelocity v of the fluid, determines the sum of the intensities due tothe volume of each phase present in the fluid. In fact, theseintensities are weighted with the polarization factors f_(α)(v) of eachphase.

In general the exponentials present in the equation (13) can beapproached by the first terms of their Taylor-series expansion whichremarkably facilitates the fit of the expression (13) to theexperimental data.

In the case where the circulating fluid is composed by two phases, byexample petroleum and water, the y-axis intercept at t=0 of the fit ofthe CPMG sequence by the expression (13) will be given by the followingexpression:

I ₀(v)=x _(p) I _(p) ⁰(v)+(1−x _(p))I _(w) ⁰(v)  (14)

where x_(p) represents the petroleum fraction, I_(p) ⁰(v) and I_(w) ⁰(v)represent, respectively, the y-intercept at t=0 for total content ofpetroleum and water, circulating at an average velocity v.

Regrouping terms in equation (14) we obtain:

$\begin{matrix}{{{I_{0}(v)} = {{x_{p}\underset{\underset{C{(v)}}{}}{\left( {{I_{p}^{0}(v)} - {I_{w}^{0}(v)}} \right)}} + \underset{\underset{D{(v)}}{}}{I_{w}^{0}(v)}}}{{I_{0}(v)} = {{x_{p}{C(v)}} + {D(v)}}}} & (15)\end{matrix}$

Taking into account the expression (15) it can be seen that, for a givenaverage velocity v of the fluid, there exists a linear relationshipbetween the y-axis intercept at t=0 of the linear fit and the fractionof petroleum x_(p) present in the fluid. So, if we carry out a fit ofthe experimental values of the y-axis intercepts at t=0 for differentknown petroleum concentrations x_(p) by the expression (19), we canobtain the coefficients C(v) and D(v). Once the coefficients C(v) andD(v) are known, it is possible to obtain the value of an unknownpetroleum concentration x_(p) by the expression:

x _(p)=(I ₀(v)−D(v))/C(v)  (16)

Determination of the fraction of components with prepolarization stageof variable length (Method 2): We define the coefficient γ_(α)(v) as thequotient of the prepolarization factors of the phase a circulating at anaverage velocity v, for two different lengths L₁ and L₂ of theprepolarization stage:

$\begin{matrix}{{\gamma_{\alpha}(v)} = {\frac{f_{\alpha}\left( {v,L_{2}} \right)}{f_{\alpha}\left( {v,L_{1}} \right)} = \frac{\left( {1 - ^{- \frac{L_{2}}{{vT}_{1}^{\alpha}}}} \right)}{\left( {1 - ^{- \frac{L_{1}}{{vT}_{1}^{\alpha}}}} \right)}}} & (17)\end{matrix}$

Applying alternatively two CPMG sequences, one with a prepolarizationstage of length L₁ and afterwards another one with a length L₂, weessentially have that the y-axis intercepts at t=0 of the linear fits ofeach pulse sequence satisfy the following equations:

I ₀₁(v)=I _(p1)(v)+I _(w1)(v)  (18)

I ₀₂(v)=I _(p2)(v)+I _(w2)(v)  (19)

In equation (19), I_(p1)(v), I_(w1)(v), I_(p2)(v) and I_(w2)(v)represent the contribution to the y-intercepts at t=0 of the linear fitsof the CPMG signals corresponding to the contents of petroleum and waterwhen a prepolarization stage of length L₁ y L₂ is used, respectively,for a fluid circulating at an average velocity v.

In addition, the following relationships must be satisfied:

$\begin{matrix}{\frac{I_{p\; 2}}{I_{p\; 1}} = {\gamma_{p}(v)}} & (20) \\{\frac{I_{w\; 2}}{I_{w\; 1}} = {\gamma_{w}(v)}} & (21)\end{matrix}$

Using the equations (18) and (19) and the equations (20) and (21) weeasily find that the values of the unknown I_(pt) e I_(w1) valuescorrespond to:

$\begin{matrix}{{I_{p\; 1}(v)} = \frac{{{\gamma_{w}(v)}{I_{01}(v)}} - {I_{02}(v)}}{\left( {{\gamma_{w}(v)} - {\gamma_{p}(v)}} \right)}} & (22) \\{{I_{w\; 1}(v)} = \frac{{{\gamma_{p}(v)}{I_{01}(v)}} - {I_{02}(v)}}{\left( {{\gamma_{p}(v)} - {\gamma_{w}(v)}} \right)}} & (23)\end{matrix}$

Knowing the calibration values I⁰ _(p1)(v) and I⁰ _(w1)(v), when eachone of the phases separately is circulating at an average velocity v ina full conduit, we can calculate the absolute petroleum and formationwater cuts circulating through the pipe, in the presence or, as well, inthe absence of gas, air or a mixture of both.

What precedes, is referred to the method of measurement of the averagevelocity in a volume element of the production conduit, previouslyselected by the introduction of a magnetic field gradient which can bepermanent or pulsed, and the selective irradiation of a volume insidethe space bounded by the excitation antenna by a pulse of a narrowbandwidth—or “soft” pulse. In order to determine the flow regime in theproduction line, it is necessary to measure the velocity profile of thecomponents of the fluid in the plane perpendicular to the direction ofthe flow, sweeping all the volume elements and applying a suitablemethod for data reduction to construct the image of the velocity profileon that plane. In the next, the measurement procedure is explained.

FIG. 9 shows the steps involved in the measurement sequence. The flow350 enters into the analyzer set following a longitudinal direction andbecomes polarized by the action of the second prepolarization field.Then, a selective excitation is performed to determine the velocityprofile of each component of the fluid. The selective excitation isperformed by the magnetic field gradient 300 along the plane which istransverse to the flow, and the excitation of the corresponding MRsignals is accomplished by a first pulse frequency selective 310. Themagnetic field gradient 300 must be confined essentially to the volumeoccupied by the radiofrequency antenna of the MR sensor and can belinear or radial with respect to that plane. The selectiveradiofrequency pulses of variable frequency 310, are modulated inamplitude and/or phase, preferably by a “sinc” function. After theapplication of the CPMG measurement sequence with space selection of theexcitation 320, the velocity of the component of the fluid into theplane or slice 330, which has been excited by the first selective pulse310 to a certain resonance frequency, is obtained. In the next step, itis selected another volume element 320 to be excited by means the pulse310 keeping fixed the magnetic field gradient 300. The average velocityof the fluid inside this volume element is measured applying thepreferred CPMG pulse sequence 320. These steps are repeated a suitablenumber of times for a fixed effective length of the prepolarizationfield 360. In the next step, the effective prepolarization length 360 ischanged to the new value 560, and the previous procedure is repeated.All the procedure just described must be repeated the number of timesnecessary in order to obtain the data matrix needed to establish: i) Thevelocity profile projected along the magnetic field gradient axis 300,and ii) The relative fractions of each component of the fluid. Finally,in order to obtain a bi-dimensional velocity profile, the magnetic fieldgradient 300 is turned by a suitable angle, preferably turning all theanalyzer set 700 and the previous sequences must be repeated. To obtainthe bi-dimensional image of the velocity profile, is preferably adoptedfor this invention the simplest method known as“projection-reconstruction”; this method is described in: P. T.Callaghan, “Principles of Nuclear Magnetic Resonance Microscopy”,Clarendon Press, Oxford (1991).

A preferred embodiment refers to the determination of the velocityprofile in a tube of cylindrical section. The magnetic field in the MRsensor module must have a radial field gradient into the transverseplane to the direction of the flow. The selective excitation willprovide, for each preferred CPMG sequence of measurements, the averagevelocity in a transverse ring, whose radius will depend on the value ofthe magnetic field and the resonance frequency of excitation. Thisprocedure is repeated for different values of the excitation frequency.

The procedure of measurement and analysis by magnetic resonance isperformed by a conventional device: a transmitter Tx, with a pluralityof filters, impedance matching and “blanking” circuits and a circuit tocontrol the quality factor Q, a receiver Rx, with a plurality ofreception channels and a plurality of filters, impedance matching and“blanking” circuits, a plurality of transmission/reception switches, aradiofrequency signals generator, preferably synthesized, a pulsegenerator, which provides a plurality of synchronized pulses of command,a device for digitizing the MR signal with a plurality of channels ofreception/digitalization and a plurality of interfaces of connection tothe control computer, to the measurement and analysis MR module and tothe prepolarization module of variable length. A medium skilledtechnician in this matter can design and implement a measurement andanalysis device of magnetic resonance for the present flow regimeanalyzer.

Where it is declared or described that a measurement pulse sequence is aCPMG it must be understood that this sequence is the preferred one forthis invention. However, an alternative embodiment for this inventioncan be declared or described as any sequence for measuring the velocityprofile of each component of the fluid which encodes or select anexcitation volume element of certain antenna and, after that, anothertemporal pulse series of reading for measuring the time of permanenceinside the MR sensor of the selected volume element be applied. Inaddition, it is declared that the pulse sequences known as SSFP and/orCPMG (described as an example in D. J. Pusiol, “Device and method forreal time direct measurement of the proportion and flow-rate of amulti-component fluid”, US 2006/0020403 A1; A. Abragam”, “The Principlesof Nuclear Magnetism”, Oxford University Press, 1973; y Bageira deVazconcelos Azeredo R., Colnago A. L. and Engelsberg M., AnalyticalChemistry, 72, 2401, 2000)), are another preferably measurementsequences.

Where it is declared or described that an apparatus of this inventionincludes, contains, posses, is compound or is constituted by certaincomponents, it must be understood, except when this declaration ordescription expresses the contrary, that one or more of that or thoseexplicitly described components can be present in the apparatus.Nevertheless, in an alternative embodiment, the apparatus of thisinvention can essentially be declared or described as consisting ofcertain components, in which the components of this embodiment couldmaterially alter the operation principle or the differentiatingcharacteristics of the apparatus, could not be present in thedeclaration or the description of this alternative embodiment. Inanother alternative embodiment, the apparatus of this invention can bedeclared or described as consisting of certain components, in whichother components of the embodiment could not be declared or described.

Where the article “a” is used in a declaration of or in a description ofthe presence of a component in the apparatus of this invention, it mustbe understood, unless this declaration or description expressesexplicitly the contrary, that the use of the indefinite article does notlimit the presence of the component in the apparatus to one in number.

1. A magnetic resonance based apparatus to analyze and to measure theuni and bi-directional flow regime in a transport or production pipe ofcomplex fluids, in real time and flow, said apparatus including ameasurement and analysis magnetic resonance module, two prepolarizationmagnetic modules of variable effective length and a control, acquisitionand data transfer computer, all mutually interconnected.
 2. An apparatusof magnetic resonance to analyze and to measure the uni andbi-directional flow regime in a transport or production pipe ofmultiphase fluids in real time and volumes, according to claim 1,wherein said measurement and analysis magnetic resonance moduleincludes: a. device that generates a homogenous constant magnetic field,thermally compensated and with addition of a magnetic field gradient,for spatial location. b. an exciting/receiving radio frequency antenna,c. a plurality of electronic modules for magnetic resonance.
 3. Amagnetic resonance apparatus to analyze and to measure the uni andbi-directional flow regime in a transport or production pipe ofmultiphase fluids in real time and volumes according to claim 1, whereinsaid prepolarization magnetic modules of variable effective lengthcomprise successive cylindrical segments of a discrete Halbach type,that can be interconnected to select the fluid components withsuccessive minor values of the spin-lattice relaxation time.
 4. Ameasurement and analysis magnetic resonance module according to claim 2,wherein said device generates an homogenous constant magnetic fieldthermally compensated with addition of space location gradient ofdiscrete Halbach type, and said device comprising a plurality ofmagnetic bars of different sizes and materials, which comprise aplurality of discs of adequate thicknesses, that are piled up andinterconnected, acting to: a. generate the main magnetic field ofmagnetic resonance, b. homogenize said main magnetic field, c.compensate the thermal variations of said main magnetic field, d.generate a gradient of reorientable local field with respect to thetransversal symmetry of the flow.
 5. A measurement and analysis magneticresonance module according to claim 2, wherein said receiving/excitingradio frequency antenna is preferably of a solenoid type when this mainmagnetic field is perpendicular to the flow direction.
 6. A measurementand analysis magnetic resonance module according to claim 2, whereinsaid receiving/exciting antenna of radio frequency is preferably of thebird cage type when this main magnetic field is longitudinal respect tothe flow direction.
 7. A measurement and analysis magnetic resonancemodule according to claim 2, wherein said plurality of electronicmodules for magnetic resonance includes: a. a transmitter T_(x), with aplurality of filter, adaptation of impedance and “blanking” circuits, b.a control circuit of the merit factor Q, c. a receiving R_(x), with aplurality of reception channels and a plurality of filter, adaptation ofimpedance and “blanking” circuits, d. a plurality of sending/receivingkeys, e. a generator of preferably synthesized radio frequency signal;f. a pulse generator, which provides a plurality of synchronized commandpulses; g. a digitalization device of the magnetic resonance signal witha plurality of reception/digitalization channels; h. a plurality ofconnection interfaces of said control computer to said measurement andanalysis module of magnetic resonance and said magnetic prepolarizationmodule of variable effective length.
 8. A measurement and analysismagnetic resonance module according to claim 2, wherein said pluralityof electronic, digital radio frequency modules for magnetic resonancecan be located to operate remotely with respect to the flow line of themultiphase fluid.
 9. A magnetic prepolarization module of variableeffective length according to claim 3, wherein said module includes aplurality of magnetic disks of discrete Halbach, which can be reorientedwith respect to each other to obtain the effect of variable effectivelength, with respect to the spin-lattice relaxation time of each phaseof the multiphase fluid.
 10. A magnetic prepolarization module ofvariable effective length according to claim 3, further comprising aplurality of discreet Halbach magnetic disks, which can be moved away orapproached to each other in the flow direction to obtain the variableeffective length effect, with respect to the spin-lattice relaxationtime of each phase of the multiphase fluid.
 11. A device that generatesa homogenous and thermally compensated constant magnetic field withaddition of space location gradient of discrete Halbach type accordingto claim 4, wherein said plurality of magnetic bars of different sizesand materials can be arranged to generate a perpendicular magnetic fieldto the flow direction.
 12. A device that generates a homogenous constantand thermally compensated magnetic field with addition of space locationgradient of discrete Halbach type according to claim 4, wherein saidplurality of magnetic bars of different sizes and materials can bearranged to generate a magnetic field along the flow direction.
 13. Adevice that generates an homogenous constant and thermally compensatedmagnetic field with addition of space location gradient of discreteHalbach type according to claim 4, wherein said plurality of magneticbars of different sizes and materials can be arranged in two concentriccylinders, one to generate the homogenous and thermally compensated mainmagnetic field of magnetic resonance, and the other to generate themagnetic field gradient, both cylinders rotatable throughout theparallel axis to the fluid vein.
 14. A device that generates anhomogenous constant and thermally compensated magnetic field withaddition of space location gradient of discrete Halbach type accordingto claim 4, wherein said homogenization of said principal field isrealized by a set of magnetic bars located in the diagonals with respectto the direction of said principal field, to compensate the quadrupolemoment, and another set of said magnetic bars in the directrix betweensaid bars generating said principal field to compensate multipolarcontributions of major importance.
 15. A magnetic resonance based methodto analyze and to measure the uni and bi-directional flow regime in atransport or production pipe of complex fluids, in real time and flow,wherein said method includes the following steps: a. the prepolarizationmodule is set in a way to obtain a zero effective length; b. of theprepolarization magnetic field, said module being located on the rightend of the measurement and analysis magnetic resonance module; c. theprepolarization module is set in a way to obtain a first effectivelength of the prepolarization magnetic field, said module being locatedon the left end of the measurement and analysis magnetic resonancemodule; d. a first magnetic field gradient orientation is selected inthe measurement and analysis magnetic resonance module; e. the selectivemeasurement of the velocity profile and the relative proportions of theselected components of said fluid is started, preferably by theapplication of the volume selective CPMG pulses sequence, according tothe following steps: i. a resonance frequency for a certain volumeelement inside the excitation and detection coil is selected by saidcontrol computer and said synthesized radio frequency generator, ii. thenuclei pertaining to a volume element are excited with a first softpulse of 90°, by said pulses generator and said radio frequencytransmitter, iii. the slope and the y-axis interception at t=0 of thespin echoes decay linear fit of the selected fluid phase and volume areread, which are obtained by application of a plurality of refocalizationand reading pulses, by said radio frequency receiver, said magneticresonance signal device, said plurality of connection interfaces andsaid control computer, iv. the slope and the y-axis interception at t=0of the spin echoes decay linear fit are stored in said control computer,v. another portion of the inner volume of the excitation coil isselected, by a suitable modification of the excitation frequency of saidfirst pulse, by said control computer and said synthesized radiofrequency generator, vi. steps i) to vi) are repeated, until exhaustingall the selected local volumes whose speed are to be measured; f.another gradient direction of the magnetic field is selected; g. steps dand e are repeated till exhausting all the non-redundant directions ofthe magnetic field gradient generated in the module; h. another lengthof the prepolarization field is selected; i. steps c to g are repeateduntil exhausting the information necessary to measure the relativeproportions of the different phases that compose the fluid and that flowfrom left to right; j. by said control computer relative proportions arecalculated for each selected volume elements, constituting amultidimensional matrix, with information of the average velocity andthe relative proportions for each selected volumes; k. by said controlcomputer, the bi-dimensional image of the velocity profile for eachcomponent of the fluid moving from the left to the right is obtained; l.to analyze the right to left flow the procedure above described is nowrepeated exchanging respectively left by right words.
 16. A velocityprofile and relative proportions of the selected components of saidfluid measurement method according to claim 15, wherein said method isperformed by a first pulsed excitation sequence to select the internalvolume element in the excitation/detection antenna and afterwards, asecond reading pulsed sequence is applied.
 17. A velocity profile andrelative proportions of the selected components of said fluidmeasurement method according to claim 15, wherein said sequence ofexcitation/detection pulses is of SSFP type.
 18. A velocity profile andrelative proportions of the selected components of said fluidmeasurement method according to claim 15, wherein said sequence ofexcitation/detection pulses is of CWFP type.
 19. A magnetic resonancebased apparatus to analyze and to measure the uni and bi-directionalflow regime in a transport or production pipe of complex fluids, in realtime and flow, wherein said apparatus includes two measurement andanalysis magnetic resonance modules, a magnetic prepolarization moduleof variable effective length and a control, data acquisition andtransfer computer, all mutually interconnected.
 20. An apparatus tomeasure and analyze the uni and bi-directional flow regime in atransport or production pipe of multiphase fluids in real time andflow-rates, wherein said apparatus can be used in both surface and downhole, as well as offshore or in the continent.